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Abstract: We study the possibility of a light Dark Matter (DM) within a constrained Min- 
imal Supersymmetric Standard Model (cMSSM) framework augmented by a SM singlet- 
pair sector to account for the non-zero neutrino masses by inverse seesaw mechanism. 
Working within a 'hybrid' scenario with the MSSM sector fixed at high scale and the sin- 
glet neutrino sector at low scale, we find that, contrary to the case of the usual cMSSM 
where the neutralino DM cannot be very light, we can have a light sneutrino DM with mass 
below 100 GeV satisfying all the current experimental constraints from cosmology, collider 
as well as low-energy experiments. We also note that the supersymmetric inverse seesaw 
mechanism with sneutrino as the lightest supersymmetric partner can have enhanced same- 
sign dilepton final states with large I^t coming from the gluino- and squark-pair as well as 
the squark-gluino associated productions and their cascade decay through charginos. We 
present a collider study for the same-sign dilepton+jets+^r signal in this scenario and 
propose some distinctions with the usual cMSSM. We also comment on the implications of 
such a light DM scenario on the invisible decay width of an 125 GeV Higgs boson. 
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1. Introduction 

While the existence of Dark Matter (DM) in our universe is by now well-established from 
various astrophysical and cosmological observations its nature and properties are still 
unknown. Many experimental efforts are under way [|], || to identify the DM candidate 
in various direct detection experiments through its scattering off different nuclei as well as 
from indirect detection through its annihilation products. Since no particle in the Standard 
Model (SM) can qualify as the DM candidate which is required to be dominantly 'cold' Q, 
these studies will also be sensitive probes of physics beyond the SM and supplement the 
new physics search at the Large Hadron Collider (LHC). 

Many extensions of the SM indeed require the introduction of new particles, some of 
which could be DM candidates. Perhaps the most popular candidate for particle DM is 
the lightest supersymmetric particle (LSP) in imparity conserving supersymmetric (SUSY) 
models Q. In the Minimal Supersymmetric extension of the SM (MSSM), the lightest 
neutralino is the usual DM candidate, as the other viable candidate, namely the scalar 
superpartner of the left-handed (LH) neutrino, is strongly disfavored by a combination of 
relic density, direct detection and invisible decay width of the SM Z-boson constraints 0]. 
More specifically, the unsuppressed coupling of the LH sneutrino to the SM Z-boson leads 
to a large cross-section for elastic scattering with the target nuclei in direct detection 
devices ||; such cross-sections are already ruled out experimentally over almost the entire 
viable mass range By the same argument, their large Z-coupling leads to far too rapid 
annihilation for the LH sneutrinos and too small a relic density compared to the Wilkinson 
Microwave Anisotropy Probe (WMAP)-measured value Q. One could make them very 
light (of order GeV) || in order to suppress the annihilation rate as well as to evade the 
direct detection bounds due to the limited sensitivity of the experiments at low masses. 
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However, a very light sneutrino is excluded by the measurement of the Z-boson invisible 
decay width at the Large Electron-Positron (LEP) collider [[H]]. 

On the other hand, the recent data from three direct detection experiments, namely 
DArk MAtter (DAMA) Q, Coherent Germanium Neutrino Technology (CoGeNT) @ 
and Cryogenic Rare Event Search with Superconducting Thermometers (CRESST) p3| , 
have suggested the hints of a light DM with mass in the 10 — 100 GeV range and cross 
section in the range 10 -3 — 10 -6 pb for elastic scattering off nucleons. Although there is 
no unanimity among these results and several other direct detection experiments j?], 14, 



15, 16, 17, 18, [L9] do not see any such positive hints of a particle DM, it has certainly 



generated considerable curiosity in a light DM scenario which can explain some/all of the 
hints of positive detection while being consistent with the null results from other direct 
detection experiments [2C]. Therefore, it may not be premature to examine some beyond 
SM scenarios accommodating a light DM candidate in case any of these positive hints are 
confirmed in near future. 

If we assume gaugino mass unification in the MSSM, the LEP collider searches of 
SUSY put a lower bound on the lightest neutralino mass of around 50 GeV [21], and the 
recent LHC data push this bound to more than about 200 GeV M, M. Even if we do 
not assume gaugino mass unification, one could derive a lower limit on the neutralino LSP 
mass of ~ 20 GeV [24] just requiring the observed DM relic density, together with the 
LEP constraints on chargino and slepton masses 1 . Therefore, if the DM indeed turns out 
to be very light as suggested by some of the recent experiments [11, 12, 13 1, we need to 
go beyond the universal scenario. Since the MSSM anyway cannot be a complete theory 
and needs to be extended to accommodate the observed small neutrino masses |26[| , it 
would be interesting to see if these extensions can also provide a viable light DM candidate 
while satisfying both the collider and relic density constraints as well as other low-energy 
constraints in the leptonic sector. 

A simple way to understand the smallness of neutrino masses is by the seesaw mecha- 
nism ||27| . The canonical type-I seesaw requires the addition of one set of heavy SM singlet 
Majorana fermions to the particle content of the SM, and the smallness of the light neu- 
trino masses are due to the heaviness of the Majorana neutrinos. In its supersymmetric 
version, the lightest superpartner of the singlet neutrinos with a small admixture of the 
left-sneutrino can be the LSP with some fine adjustment of parameters or in an extension 
of the MSSM gauge group (see e.g., Ref. ]^8| and references therein). In such a case, it can 
be a viable light DM candidate. 

In this paper, we focus on the supersymmetric version of a different seesaw scenario, 
namely, the inverse seesaw ]^9| where one adds two sets of SM singlet fermions, one Dirac N 
and one Majorana S per family, to generate the small neutrino masses. These models have 
three lepton-number carrying electrically-neutral fermions per family, namely (vl, N c , S). 
If a linear combination of the super-partners of these fields turns out to be the LSP, then it 
could be a scalar DM candidate. Current literature on the subject discusses several varia- 
tions of such models, e.g., (i) within the framework of MSSM |p0| , (ii) in Next-to-Minimal 

For a recent discussion on how to avoid this lower bound in a non-universal MSSM parameter space, 
see e.g., Ref. p5| and references therein. 
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Supersymmetric SM (NMSSM) Q, (iii) with extended gauge symmetry SU(2)l x U(1)y x 
?7(1)b-l (3^] so that the seesaw mass matrix arises from a B — L gauge symmetry |33|] 2 , 



(iv) with Supersymmetric Left-Right gauge group SU{2)l x SU(2)r x U(1)b-l [ 36 1 where 
the inverse seesaw texture arises naturally even at TeV scale while preserving the gauge 
coupling unification [37]. 



In the present work, we take a hybrid approach similar to that in Ref. |3Cf| , i.e., a low- 
energy input for the SU (2) L-singlet neutrino sector and for the lepton-number violating 
soft SUSY-breaking sector while a top-down approach for the MSSM particle spectrum, 
without necessarily imposing any features of a specific Grand Unified Theory (GUT)-based 
model. Our goal is to examine if such a minimal Supergravity (mSUGRA) scenario with 
inverse seesaw can give a light DM candidate satisfying all the existing cosmological, collider 
as well as low-energy constraints. If so, we ask ourselves what the collider signals for such a 
scenario are and how to distinguish it from a typical constrained MSSM (cMSSM) scenario 
for a similar squark-gluino spectrum at the LHC. In particular, since the sneutrino mass is 
not directly related to the gaugino masses, and there are additional unknown parameters in 
the sneutrino mass matrix, we expect the lightest sneutrino to be allowed to have masses 
in the few GeV range without being in conflict with the collider bounds on gluino and 
chargino masses. Then it remains to be checked whether this lightest sneutrino eigenstate 
has the right admixture of left- and singlet-sneutrino flavors to reproduce the observed relic 
density while satisfying the constraints from direct and indirect detection experiments as 
well as from other low-energy sectors. In fact, we find that in contrast with the usual 
cMSSM scenario p5| 1, we can have the sneutrino LSP mass in the few GeV range while 
being consistent with the SUSY search limits; however, the relic density constraint, among 
others, requires the lightest sneutrino mass to be more than ~ 50 GeV. Though this is not 



consistent with the CoGeNT-preferred range of ~ 10 GeV 12] for the DM mass, it is within 



the 2ci-preferred range of CRESST-II [ 13 1 and also close to one of the DAMA-preferred mass 



range j[l]]. Moreover, the benchmark points we find around 50 GeV sneutrino DM mass 



are all consistent with the recent hints of the lightest Higgs mass around 125 GeV 38, 39] 



which is very difficult to accommodate for a light neutralino DM in the usual cMSSM 



scenario |23|, |40 |. 



This paper is organized as follows: in Section 2, we briefly discuss the SUSY Inverse 
Seesaw Model (SISM) parameters and set up our notation; in Section 3, we present some 
benchmark points for the sneutrino DM in SISM satisfying all the existing experimental 
constraints; in Section 4, we discuss the collider signatures of a sneutrino LSP through cas- 
cade decays of squarks and gluinos and identify the same-sign dilepton+jets+large $t sig- 
nal; in Section 5, we present a detailed collider simulation of this signal for our benchmark 
points, along with the relevant SM background, and propose some distinctions with the 
usual cMSSM case having similar squark-gluino spectrum; and in Section 6, we summarize 
our results. 



2 Another class of models [Q uses global B — L symmetry to restrict the inverse seesaw matrix to the 
desired form. Yet another recent work [^] uses a gauged [/(l)-extension with only renormalizable operators, 
but adding pairs of fermions with fractional lepton number. 
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2. The Model 



In order to explain the non-zero neutrino masses by an inverse seesaw mechanism |29|| , the 
MSSM field content is supplemented by three pairs of SM-singlet superfields, (Dirac) iV, 
and (Majorana) Si (i = 1,2,3 for three generations) 3 . The superpotential is given by 

Wsism = Wmssm + eaby'jL^NjHt + M R%] N.Sj + SiSj (2.1) 

where the /is-term is the only lepton-number breaking term in the superpotential. The 
corresponding soft SUSY-breaking Lagrangian is 



/'soft /'soft 

^SISM — ''-MSSM 



+ rriq S^S 



e ab ^L1N 3 H h u + B^NiSj + B%S % S + h.c. 



(2.2) 

The tree-level 9x9 neutrino mass matrix in the basis {vli N c , S} is given by 

/ Md 

Ml M R I (2.:!) 
V Ml us 



Mi 



where Mb = v u y v is the Dirac neutrino mass matrix, v u being the vacuum expectation 
value (vev) of the -ff u -superfield. The 3x3 light neutrino mass matrix in the approximation 
fis *C Md < Mr is given by 

M v = {MdM^ 1 ) fis {MoM^f = Ffi s F T (2.4) 

Assuming a TeV-scale inverse seesaw, Mr ~ O(TeV), and O(0. 1) Dirac Yukawa coupling 
(i.e., Md ~ O(10) GeV), we need the lepton-number violating mass term us ~ O(keV) for 
a sub-eV light neutrino mass, as required by the neutrino oscillation data |26[ |. Here we 
emphasize the fact that us is much smaller than the other energy scale(s) pertinent to the 
SUSY sector. This feature, namely the lepton number violation at a very low scale, is the 
quintessence of the inverse seesaw mechanism, which is integrated with the SUSY scheme 
here. The smallness of us is technically natural in the 't Hooft sense, but must have its 



origin from some other new physics, e.g., radiative corrections jl^] or extra dimensions [43]. 

We note here that the gauge symmetry SU(2)l x U(1)y allows for additional entries 
in the singlet sector, i.e., non-zero v^S- and iViV-ternis in the neutrino mass matrix given 
by Eq. fl2,3| ). However, the presence of only the NN term does not spoil the inverse 
seesaw structure at tree-level since the rank of mass matrix still remains the same. But 
the vlS term will, in general, affect the inverse seesaw formula given by Eq. (2.4), unless 
the coupling ys in the corresponding ysLH u S-term in the superpotential is ^ 10~ 12 or so. 
These issues can be naturally eliminated by extending the SM gauge group so that these 
additional terms in the superpotential are forbidden by some symmetry (see, for instance, 



Refs. [33, ^7|). However, in order to allow a direct comparison with the usual cMSSM case, 



we choose to work within the MSSM gauge group SU(2)l x U(1)y and assume ys = 0. 



3 Only one pair of S , [/(2)i-singlets is sufficient to satisfy the neutrino oscillation data EJ]. However, if 
we want to generate all the neutrino masses at the tree-level, we must have three pairs of singlets. 
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The mixing in the light neutrino sector is usually described by the unitary Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix U which diagonahzes the light neutrino mass ma- 
trix: 

U T M U U = diag(mi, m 2 ,m 3 ) (2.5) 

Since the above diagonalization of M v does not simultaneously diagonalize the other mass 
matrices Mr and /ig appearing in the full neutrino mass matrix given by Eq. ([2.3;), there 
will be, in general, additional non-unitary contributions to the light neutrino mixing matrix 
due to its mixing with the heavy neutrinos. This can be derived from the 9x9 unitary 
matrix V which diagonahzes the full neutrino mass matrix given by Eq. (|2.3|): 



VM V V T = diag(m i ,m JJ .), (i = 1, 2, 3; j = 1, 2, .., 6) (2.6) 
by decomposing it into the blocks 

v 9x9 =f"r v x6 V (2 - r) 

^6x3 -^6x6 J 

Then the upper 3x3 sub-block will represent the full (non-unitary) light neutrino mixing 
matrix. To leading order in F = MpM^ 1 , this can be expressed in terms of the PMNS 
matrix as follows: 



U~(i-~FfAu=(1-tj)U (2.8) 

where r] = \FF^ measures the non-unitarity of the light neutrino mixing matrix. 

In the corresponding scalar sector, the sneutrino mass matrix is a 9 x 9 complex, or 
18 x 18 real matrix which can be decomposed into two 9x9 block-diagonals assuming 
CP-conservation in the corresponding soft-breaking sector: 

J with 

/ m| + M D M ] D + \m\ cos 2/3 ±(v u A u - p,M D cot /3) M D M R 

±{v u A u - fiM D cot (3) m 2 N + M R M R + M D M j D ix s M R ± B Mr 

\ M D M R h s M r ±B Mr m 2 s + n 2 s + M R M R ±B fJ 




Ml = 



where the corresponding mass eigenstates are linear combinations of the three sneutrino 
flavor eigenstates: vr+j. = (vL a , N£, Sd^j = 1,2,. ..,9; a,b,d = 1,2,3). In the next sec- 
tion, we examine the SUSY parameter space in which the lightest of these mass eigenstates 
can be the LSP. 

3. Some Benchmark Points 

Our goal in this section is to find a sparticle spectrum with light sneutrino LSP in the 
cMSSM scenario with 5 parameters (mo, m 1 / 2 , tan (3, Aq, sign fj) and the additional inverse 
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Table 1: Various experimental constraints used in our analysis to find the benchmark points. 

seesaw parameters /is, Mr, Mjj, and Bm r - Once we find a light sneutrino LSP, we 
impose the relic density and direct detection constraints in order for it to qualify as a 
DM candidate. We also require all the benchmark points to satisfy various collider and 
low-energy constraints, summarized in Table [|. A few comments: 

• For the lightest Higgs mass, we use the CMS suggested value of 125.3 ± 0.4(stat) ± 
0.5(syst) @. The ATLAS suggested central value is around 126.5 GeV || with 
presumably similar experimental uncertainties; for concreteness, we just choose to 
work with the CMS value which has the errors explicitly stated. 
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For the spin-independent DM-nucleon scattering cross-section, we use the 2a upper 
limit from the latest XENON100 data @. 

There also exist strong constraints on DM annihilation cross sections from indirect 
detection searches, e.g., in gamma rays 47] and in high-energy neutrinos from 
the Sun [53|. Here we use the latest 95% confidence level upper limits obtained from 
the Fermi Large Area Telescope data [47]. 

The lepton anomalous magnetic moments as shown in Table || are defined as Aae = 
o| M — a^ xpt where ap = (g — 2)^/2. The most important one is the muon anomalous 
magnetic moment which persistently shows a 3cr discrepancy [S3] over the SM pre- 
diction and should be taken into account in any complete beyond SM scenario. For 
the electron (<? — 2), the discrepancy is quite small and is a rather loose constraint on 
the new physics parameter space. We do not consider the tau anomalous magnetic 
moments here, because its value is not known so precisely 0. 



The non-unitarity of the light neutrino mixing matrix is defined in Eq. fl2,8|) and the 
constraints on its elements shown in Table [T| are derived from a combination of the 
neutrino oscillation data, the LEP precision data from weak gauge boson decays and 
the lepton-flavor violating (LFV) decays p2[ . 



There are also strong constraints on the cMSSM parameter space from direct SUSY 
searches at the LHC [54, |55|. Therefore, we must choose the input points in the (mo, mi/2)- 



plane not already excluded by the LHC SUSY searches which for certain cases extends to 
m i/2 ~ 600 GeV and mo ~ 1 TeV (e.g., in the jets+^r channel |56|]). On the other 
hand, very large values of mo and mi/2 (larger than a few TeV) are not desirable from 
phenomenological perspective as they drive most of the sparticle masses beyond the kine- 
matic reach of the LHC. Therefore, we choose our mo values close to 1 TeV and the m 1( / 2 
values close to 600 GeV. We also choose to work with /x > case, since [i < is strongly 
disfavored by the muon anomalous magnetic moment as well as by the B — > X s ^f branching 
ratio f57f| . Similarly, large tan f3 values ^ 50 are disfavored by the recent LHCb results on 
B s — > /x + /i _ f5l"[] , and hence, we choose some intermediate values between 25 and 35 for the 
benchmark points discussed below. For the trilinear term Aq, the recent LHC discovery of 
a SM Higgs-like particle at 125 GeV |3q , |39j ] implies that we must have a large negative 
A-teim (for /x > 0) in order to have the radiative corrections account for the required 
enhancement of the lightest Higgs mass from its tree level value close to mz |]58| . 

In the neutrino sector, for simplicity, we assume the inverse seesaw parameter matrices 
Md,Mr as well as the l?-terms B fls ,BM R to be diagonal 4 . Hence, we can easily satisfy the 
LFV constraints for our benchmark points. Allowing non-zero off-diagonal entries in the 
Dirac Yukawa coupling matrix y u will induce large LFV effects, and we find that for the 
benchmark points discussed in the following section, we must have the off-diagonal entries 
less than ~ 0(0.01) in order to satisfy all the LFV decay modes listed in Table ||. Moreover, 
we assume no CP-violation in the neutrino sector, and choose all the mass matrices to be 



4 We can choose this kind of texture since we are not working within any particular GUT framework. 
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real 5 . Fixing both Mjj and Mr also fixes the lepton-number breaking Majorana mass 
matrix us by fitting to the neutrino mass and mixing parameters (assuming a particular 
mass hierarchy for the light neutrinos). Also note that since we are assuming a complete 
unification of the scalar sector, we choose m\ = m? N = m s = uiq at the high scale and 
similarly for the A-terms. 

The input parameters are chosen in such a way that all the experimental constraints 
listed in Table |l] are satisfied for all the benchmark points. Table ^ lists all the input 
parameters for three benchmark points we have chosen to work with. For the low-energy 
values of y u and Mr obtained by the renormalization group evolution of the parameters 
given in Table |2j, the observed neutrino mass and mixing parameters can be fitted using 
appropriate values for the mass matrix /is in Eq. fl2.4| ). As an example, for a normal hier- 
archy of light neutrino masses, using the latest global fit values for the neutrino oscillation 
parameters |5^] which includes the most recent #13 results from Double CHOOZ, Daya Bay 
and RENO experiments: 

Am^ = (7.62 ± 0.19) x 10~ 5 eV 2 , Am 2 31 = (2.53 ± 0.09) x 10" 3 eV 2 , 

sin 2 9 12 = 0.320 ± 0.016, sin 2 6» 23 = 0.490 ± 0.065, sin 2 6 13 = 0.026 ± 0.004, 



we obtain the /is values as shown in Table 



Input parameter 


BP1 


BP2 


BP3 


m (GeV) 


993.68 


996.84 


815.79 


m 1 / 2 (GeV) 


600 


650 


600 


A (GeV) 


-2712.11 


-2858.42 


-2442.11 


tan /3 


35 


25 


30 


Vu 

Mr (GeV) 
B w (GeV 2 ) 
B Mr (GeV 2 ) 

MS (eV) 


(0.16,0.16,0.18) 
(300,1000,1000) 

10 

10 6 

/ 2.04 8.27 -5.34 \ 
1 8.27 56.69 12.31 J 
V -5.34 12.31 79.26 / 


(0.10,0.10,0.08) 
(200,1000,1000) 

10 

10 6 

/ 2.25 12.82 -9.69 \ 
j 12.82 123.64 31.39 ] 
V -9.69 31.39 236.39 / 


(0.10,0.10,0.10) 
(610,1000,1000) 

10 

10 6 

/ 19.77 38.01 -23.11 \ 
I 38.01 123.87 25.31 
V -23.11 25.31 153.34 / 



Table 2: The input parameters for three chosen benchmark points (BP). The mSUGRA parameters 
are defined at the high scale and the singlet neutrino parameters at the low scale. We assume fi > 
throughout and the neutrino sector parameters shown here have been chosen to be diagonal, as 
discussed in the text. 

The low-energy mass spectrum for the superpartners corresponding to the three bench- 
mark points are tabulated in Table ||. For this purpose, we have used SARAH |Mj to imple- 
ment the SISM scenario, and SPheno |il| to generate the mass spectra and to evaluate some 
of the low-energy observables. The DM relic density and its scattering and annihilation 
cross sections were calculated using micrOMEGAS |62|| . 

5 The addition of one or more CP-phases in the neutrino sector will not affect the sparticle spectrum, 
and hence, is irrelevant for our main results in the subsequent sections. 
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1247.1 


1125.3 




d x 


1531.6 


1603.4 


1430.2 




^2 


1571.1 


1647.4 


1473.1 




Sx 


1531.5 


1603.3 


1430.1 




s 2 


1570.9 


1647.3 


1473.0 




&1 


1087.8 


1194.3 


1061.8 






1304.0 


1460.0 


1265.4 


Gluino 


5 


1401.4 


1505.3 


1392.6 




X? 


264.3 


286.2 


261.8 


Neutralino 


X% 


499.2 


539.8 


495.2 




X°3 


-1376.5 


-1464.4 


-1295.3 




xi 


1379.5 


1467.5 


1298.7 


Chargino 


xt 


499.4 


540.0 


495.4 




Xt 


1380.1 


1467.9 


1299.2 



Table 3: The sparticle masses (in GeV) for the chosen benchmark points. The sneutrino masses are 
shown up to four decimal places to illustrate the lifting of degeneracy between the mass eigenstate 
pairs due to the small lepton-number breaking. 

Note that the sneutrino real scalar fields , Vr. ) are split in their masses with the 
mass splitting in the range of keV-MeV within each pair which is a characteristic feature 
of the SUSY inverse seesaw mechanism [ |36[| . In later sections, we will sometimes denote 
the lightest mass eigenstate pair (uj 1 , ) simply by u\ . 

It is clear from Table || that all the benchmark points satisfy the LHC direct search 
limits on the SUSY particle masses in cMSSM. They also satisfy the other low-energy 
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experimental constraints in Table |], as shown in Table ||. Here we want to make some 
comments on these observables: 



Parameter 


BP1 


BP2 


BP3 


Tflh 


123.9 


123.8 


123.7 


a "JJJVL' 6 


0.105 




0.106 




0.119 






8.2 x 10" 


Q 

-y 


1.3 x 10 




— y 


2.4 x 10" 


Q 

— y 


(ct ait) (rm^s^^ 


5.6 x 10" 


-34 


1.5 x 10" 


-34 


5.4 x 10" 


-35 


Aa 


5.1 x 10" 


1 n 
1U 


3.3 x 10 _ 


- 1U 


5.4 x 10- 


-1U 




1.2 x 10" 


1 A 


7.7 x 10" 


-lo 


1.3 x 10- 


1 A 
-14 






-4 


1 \y 1fl" 

Z. i X 1U 


-4 


Z.D X 1U 


-4 


RR ( R _i. + 


3.6 x 10 


-4 


3.7 x 10" 


-4 


3.7 x 10" 


-4 


BR(/x -»■ 67) 


8.0 x 10" 


-21 


8.2 x 10" 


-22 


4.4 x 10" 


-21 


BR(r -> e7) 


1.0 x 10" 


-19 


1.6 x 10- 


-20 


8.1 x 10" 


-20 


BR(t -> p 7 ) 


1.4 x 10" 


-15 


2.4 x 10- 


-16 


1.3 x 10" 


-15 


BR(/i -> 3e) 


6.7 x 10" 


21 


2.2 x 10- 


-21 


6.2 x 10" 


-23 


BR(r -> 3e) 


1.7 x 10" 


19 


4.5 x 10- 


-20 


1.1 x 10" 


-21 


BR(r -> 3/it) 


3.6 x 10" 


15 


6.2 x 10" 


-16 


4.8 x 10" 


-17 


heel 


4.3 x 10" 


-3 


3.9 x 10 


-3 


4.5 x 10" 


-4 




3.4 x 10" 


-4 


1.6 x 10 


-4 


1.6 x 10" 


-4 


|*7tt| 


2.9 x 10" 


-4 


9.8 x 10 


-5 


1.5 x 10" 


-4 



Table 4: The low-energy observables for the three chosen BPs. These values are to be compared 
with the experimental values in Table 1. 



• It is well known that a 125 GeV mass for the lightest neutral Higgs boson in MSSM 
is not very natural |33|. It becomes even more difficult in cMSSM if one has to 
satisfy the other low-energy constraints and requires the neutralino LSP to have the 
observed relic density p3|, [40[|. The situation is somewhat similar in our case; however, 
since fixing the exact mass of the suspected scalar resonance at the LHC will require 
more data, we are content with values within 2 GeV of the average of the CMS and 
ATLAS central values. We believe that any tweaking of parameters to confirm the 
exact Higgs mass, when it is known with greater precision, will not affect the general 
conclusions of this paper. 

• The correct relic density is obtained near the resonant enhancement region of the 
annihilation cross-section in the Higgs-mediated s-channel process: — > ff where 
/ denotes the SM fermion (mostly b and r in our case). Therefore, all our benchmark 
points have the LSP mass close to m/j/2. This is illustrated in Figure [l] where we have 
plotted the relic density versus the sneutrino LSP mass for some typical Dirac Yukawa 
coupling values. We also find that the sneutrino LSP-nucleon spin-independent elas- 
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Figure 1: The relic density as a function of the sneutrino LSP mass for the first generation Dirac 
Yukawa coupling values of 0.10 and 0.16 while keeping all other input parameters fixed. The points 
corresponding to our benchmark points are shown by dots. The shaded region is the 3cr-range 
allowed by the WMAP data. 

tic 6 cross section is within the la upper bound of the latest XENON100 result except 
for BP1 which is slightly above their 2a limit. 

• The dominant annihilation channels for the sneutrino DM in our case have bb and rf 
final states. For small admixtures of the left sneutrino component (allowed by the 
Z-invisible decay width) in a mostly dominant singlet sneutrino LSP eigenstate, the 
thermally averaged annihilation cross section values for our benchmark points are 
well within the experimental upper limits given in Table |l|. 

• The SUSY contributions to Aa M have been calculated to the two-loop order |M| and 
it has been shown that for a range of parameter space, it is possible to accommodate 
the discrepancy. As can be seen from Tables [l] and |], we are able to explain the 
discrepancy within 3a for all our benchmark points. Also for the electron anomalous 
magnetic moment, we are consistent with the smaller discrepancy to within la. 

• The SM prediction for the branching ratio of the weak radiative S-meson decay is 
(3.15±0.23) x 1CU 4 J65| , and comparing with the most recent experimental value from 
the BaBar experiment as given in Table [fl, we see that there is very little room left 
for the SUSY contribution pq ]. However, for the choice of our mSUGRA parameters, 
the SUSY contributions are negligible and the values predicted for all our BPs are 
within 2a of the experimental value. On the other hand, for the branching ratio of the 
flavor-changing-neutral-current (FCNC) process — > fi + fi^, the SM contribution 
is small: (3.2 ± 0.2) x 10~ 9 |37j; hence, comparable SUSY contributions |6£| are 
still allowed by the latest data from the LHCb experiment. We estimated that the 
branching ratios for our benchmark points are within this allowed range. 

6 Note that the splitting between the lightest sneutrino mass eigenstates for our benchmark points is of 
order of a few MeVs which is too large for inelastic DM-nucleon scattering to occur. Hence we only consider 
the elastic scattering. 
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• For the rare LFV decays, since we are working within an mSUGRA scenario, the 
SUSY contributions are quite small |}{|. But the contributions from the leptonic 
sector could be large in seesaw models with large Yukawas |70fl , as in our case. How- 
ever, due to our choice of the diagonal textures for the Dirac Yukawas, the leptonic 
contributions also vanish altogether. Hence, we have very small LFV branching ratios 
for all the benchmark points. 

• The non-unitarity effects could also, in principle, be large in low-scale inverse seesaw 
models with large Dirac Yukawas |37], Ell ^TJ. In our case, again due to the diago- 
nal textures chosen for both Md and Mr, the non-unitarity parameter, defined by 
Eq. Q2.8D , is also a diagonal matrix. Hence we only show the values for its diagonal 
entries in Table |4|, and all our values are within the current experimental bounds. 

We also note that in our scenario, since the sneutrino LSP is sufficiently light, the 
lightest neutral Higgs boson can, in principle, decay into a pair of LSP's, thus giving rise 
to an invisible decay width of the Higgs boson. The LHC signatures of these decays are 
relatively clean, and very large branching ratios to an invisible decay channel are disfavored 
by the current LHC Higgs searches |[72|]. The branching ratio depends, among other things, 
on the neutrino Yukawa coupling y v . Recent global analyses [Q 7 have reported that the 
present LHC Higgs data can indeed accommodate an invisible branching ratio for the Higgs 
boson, although their best fit values for this do not quite agree with each other. If such a 
possibility is more precisely fixed by future data, it may lead to an estimate of the bounds 
on the neutrino Yukawa couplings in the inverse seesaw models which could be compared 
with those obtained from direct Higgs search results |75f| . 

We have thus demonstrated convincingly that (a) a hybrid scenario for the origin of 
soft SUSY-breaking masses can be used consistently with the inverse seesaw mechanism, 
(b) one can have a sneutrino LSP which is light and is still consistent with all the existing 
experimental constraints, and (c) the rest of the SUSY spectrum is phenomenologically 
viable. The next question to ask is whether there are any distinctive signatures of this 
scenario which can be seen at the LHC. We address this question in the next section. 



4. Collider Signatures 

The most copious collider signals of any SUSY scenario will come from the production of 
colored superpartners, namely squarks and gluinos, which will have cascade decays through 
charginos and neutralinos, eventually ending up in the stable LSP in imparity conserving 



SUSY models |76 ]. Unless these squarks and gluinos are too heavy to be kinematically 
accessible, they will have substantial production cross sections at a hadron collider due to 
the strong interaction. The production channels are gluino-pair production, squark-gluino 
associated production and squark-squark pair production (see Fig. ||). As the direct decay 
of the squarks and gluinos to the color- and electrically-neutral LSP are either forbidden or 
occur with only a tiny branching fraction, the dominant decay modes for the gluino always 



7 For a similar analysis with the earlier LEP/Tevatron/XENON/WMAP data and a 50-60 GeV scalar 
DM scenario (as in our case), see Ref. 
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involve quarks (and hence multiple jets in the final states). The gluino can have either 
the two-body decay via g — > qq, if kinematically allowed, or the three-body decay modes 
g — > qq'xf, QQX^j with virtual squarks. Similarly, the squarks decay to two-body modes 
q — > qg, if kinematically allowed, or qi — > q'xf, QX^ while <jr —> qx°j only, since right- 
handed squarks do not couple to charginos in the MSSM. If the squarks are degenerate, and 
the Yukawa coupling effects negligible, the three-body decays to the wino-like charginos 
and neutralinos usually have larger branching fractions due to their larger gauge couplings. 
If \n\ < M%, gluinos and squarks may thus decay most of the time to the heavier charginos 
and neutralinos, resulting in lengthier cascade decay chains than those shown in Figure 




Figure 2: The Feynman diagrams for the gluino pair-production and its cascade decays to give 
the same-sign dilepton+jets+^T signal at the LHC. Similar diagrams exist for squark-gluino and 
squark-squark production and decay which we have not shown here. 

The LHC signals of our scenario can differ considerably from those of the usual cMSSM 
situation. To understand this, let us look at the branching ratios of the two-body decays 
of the lighter chargino (xf) which are listed in Table [| for all of our benchmark points. 
While the conventionally expected decay xf ~^ ^ + Xi is there, it is dependent on the 
W3 components of Xi as well as the Higgsino components of both xf an d Xi- O n the 
other hand, here we have another often dominant channel, namely xf ~ * (where v\ 
is the sneutrino LSP and I = e, fi) triggered by the large mixing in the sneutrino sector 8 . 
Consequently, the leptonic branching ratio of the xf is remarkably enhanced 9 . Thus the 
SUSY cascades lead to a highly boosted rate of dileptons, of which the same-sign dileptons 
(SSD) are more spectacular being relatively background-free. The scenario outlined by 

s £ is mostly electron for our choice of benchmark points, though cases with muons do not make any 
difference in our analysis. Since the tau-lepton detection efficiency is not as good as for the electron and 
muon, we will not analyze the tau-lepton final states. 

9 Note that all charged sleptons are heavier than the lighter chargino in our case. 
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us will therefore exhibit a rise in the SSD rate with respect to that of purely jets+^T 
events, as compared to a cMSSM spectrum of comparable heaviness. We also expect the 
Jfjj> distribution to be different for a sneutrino LSP case than the neutralino LSP case, as 
noted earlier in Ref. [77]. In particular, the ^-^distribution is expected to be much harder 
in our SISM scenario compared to the cMSSM scenario. 



Chargino (Xi ) decay 


BP1 


BP2 


BP3 




0.23 
0.77 


0.45 
0.55 


0.31 
0.69 



Table 5: The lighter chargino decay branching ratios for our benchmark points in SISM. On the 
other hand, in the cMSSM case with neutralino LSP, the branching ratio is close to 100% for the 
decay W + x?. 

To illustrate this SSD-enhancement effect in our case compared to the pure cMSSM 
scenario, we construct a ratio as follows: 



a(£ ± e ± + > 2j + ftr) 
a(0£+ > 3j + fa) 



(4.1) 



which is expected to be larger in our case, and as shown in the next section, could be used 
to distinguish our SISM scenario with sneutrino LSP from the usual cMSSM scenario with 
neutralino LSP. 

Here we want to emphasize that the SSD signal in inverse seesaw is purely supersym- 
metric in nature. In other words, if one leaves aside the SUSY processes, the SSD signal 
is suppressed due to the small lepton-number violation and pseudo-Dirac nature of the 
singlet neutrinos. In that case, however, one can look for the tri-lepton signals with ^-dor 
its LHC discovery potential [78, f7{|. Note that one can also investigate the SUSY inverse 



seesaw in the tri-lepton channel in which case novel correlations of the tridepton signal 



with the neutrino mixing angles can be searched for [80]. 



5. Event Generation, Background Simulation and Results 

In this section, we give a detailed description of the SSD+jets+^T signal in our SISM 
case with light sneutrino LSP and a comparison of the signal strength with a canonical 
cMSSM scenario with neutralino LSP having similar squark-gluino spectrum for a possible 
distinction of the two cases at the yfs = 14 TeV LHC. The SUSY spectrum and the various 
decay branching fractions were calculated using SPheno pi] . The input files are then fed 
to PYTHIA (version 6.409) ]^l|] for event generation. The initial and final state quark 
and gluon radiation, multiple interactions, decay, hadronization, fragmentation and jet 
formation are implemented following the standard procedures in PYTHIA. The factorization 
and renormalization scales are set at (i.e hr = [if = V? ), where is the parton 
level center of mass energy. We have used the leading order CTEQ5L parton distribution 
functions [82] for the colliding protons. The jets are constructed using the cone algorithm 
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in PYCELL; only those jets are constructed which have pt > 20 GeV and \rj\ < 2.5. To 
simulate the detector effects, we have taken into account the smearing of jet energies by a 
Gaussian probability density function of width a(E)/Ej = (0.6/ Ej [GeV]) + 0.03 where 
Ej is the unsmeared jet energy [33|. 



In order to find the same-sign di-leptons+n jets+^T (with n > 2) final states, we 
impose the following selection criteria: 

• p T > 10 GeV and \rf \ < 2.4 for both the leptons. For the same-flavor dilepton final 
states, we raise it to p T > 15 GeV. 



• Lepton-lepton separation ARu > 0.2, where AR = \J (At/) 2 + (Acfr) 2 . 

• Lepton-jet separation AR^ > 0.4. 

• The sum of Ex deposits of the hadrons which fall within a cone of AR < 0.2 around 
a lepton, must be less than 0.2 p T . 

• Jet-jet separation ARjj > 0.4. 

Since our goal is to distinguish the SUSY inverse seesaw scenario from the conventional 
cMSSM case, we need to consider similar squark-gluino spectrum for both the cases. In 
order to do so, we generated similar benchmark points for the cMSSM case using the 
same mSUGRA input parameters given in Table || and also checked that the effective mass 
distributions, defined as the scalar sums of the lepton, jet and missing transverse energy: 



M eff = ^|^ T |+^|p J T |+^ T , 



(5.1) 



are similar for both the scenarios, as shown in Fig. |3] for all the benchmark points. 

Now in order to distinguish the two scenarios, we compute the ratio r defined in 
Eq. fl4.1|) for both SISM and cMSSM cases which are tabulated in Table |6[ We find that 
the value of r in the SISM case is roughly 4 — 5 times higher than the cMSSM case for 
all the benchmark points. Apart from this clear distinction, we also expect more in 



r 


BP1 


BP2 


BP3 


SISM 
cMSSM 


0.19 
0.04 


0.15 
0.03 


0.11 
0.03 



Table 6: The r values for all the benchmark points in both SISM and cMSSM cases. 

the chargino decay in case of the SISM, as mentioned earlier. In order to illustrate this, 
we need to analyze the fx distribution for the SSD+> 2j + fx signal for both SISM and 
cMSSM cases. We also need to analyze the SM backgrounds in detail, as follows. 

The dominant SM background for the SSD events come from tt, Wti, WWW, WWnj, WZnj, 
ZZ, Wbb, Zbb final states at the LHC |78|| . All the SM backgrounds except the tt were gen- 
erated at the parton level using ALPGEN (version 2.14) (84|] with default factorization and 
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Figure 3: The effective mass distribution of the final states as defined in Eq. (5.1) for the SISM 
and cMSSM scenarios to illustrate that the squark-gluino spectrum considered in both cases are 
similar. 

renormalization scales, and then fed to PYTHIA for showering, hadronization, fragmenta- 
tion, decay, etc. The ti background was directly generated and analyzed in PYTHIA. The 
number of events obtained after the selection criteria for yfs = 14 TeV LHC and normal- 
ized to 30 fb _1 luminosity are shown in Table [5]. Note that at this stage, some of the SM 
backgrounds are much larger than the SSD signal, and we need to devise further cuts to 
reduce the background without affecting the signal much. As shown in Table [?], we found 
two relevant cuts, namely (a) > 50 GeV for all jets and pip > 100 GeV for the leading 
jet, and (b) $t> 300 GeV which reduce the SM background significantly. 

The $t distributions for both SISM and cMSSM cases are shown in Figure |] for all 
the benchmark points. It is clear that the SISM case has a much harder tail compared 
to the cMSSM case which can be used as a distinguishing feature. The combined SM 
background is also shown (in shades) which falls rapidly for Ifi? > 300 GeV. This justifies 
our fix cut selection in Table [7|. 

6. Summary and Conclusions 

We have considered a SUSY scenario that accommodates the inverse seesaw mechanism of 
neutrino mass generation via a small AL = 2 mass term (of the order ~ keV). Two sets 
of S'C/(2)i-singlet neutrino superfields have been introduced for this purpose. We show 
that this model can not only account for the neutrino masses and mixing, but also leads 
to an LSP dominated by right chiral sneutrino states. For phenomenologicaily consistent 
input parameters, taken as a hybrid of the top-down and bottom-up choices, the sneutrino 



- 16 - 



Channel 


After basic selection criteria 


After jet-pT cut 


After $ T 


cut 




Mi 


efi 


ee 


nn 


e// 


ee 


Mi 


e/./. 


ee 


BP1 


33.24 


125.18 


144.01 


30.73 


112.66 


127.45 


24.30 


90.13 


114.69 


BP2 


39.95 


32.44 


97.26 


34.38 


26.86 


84.34 


28.54 


23.35 


64.87 


BP3 


35.94 


88.94 


102.84 


34.15 


80.05 


91.49 


32.44 


78.48 


86.76 


WWW 


16.86 


12.36 


29.64 


3.18 


2.49 


6.00 


0.39 


0.24 


0.24 


WWjj 


140.01 


75.39 


193.86 


75.39 


43.08 


96.93 


0.00 


0.00 


0.00 


wz 


84.60 


16.92 


186.06 


33.84 


0.00 


51.00 


0.00 


0.00 


0.00 


zz 


0.33 


0.33 


0.66 


0.000 


0.000 


0.03 


0.00 


0.00 


0.00 


Wbb 


29.25 


5.85 


29.25 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Wti 


81.33 


66.84 


147.54 


38.70 


31.89 


69.75 


1.83 


1.59 


3.18 


it 


2109.00 


754.80 


2331.00 


710.4 


222.00 


466.2 


0.00 


0.00 


0.00 


Zbb 


0.00 


6.99 


19.38 


0.00 


0.00 


1.62 


0.000 


0.000 


0.000 



Table 7: The number of events for 30 ftr 1 luminosity at y/s = 14 TeV LHC for the SSD+?ij + $ T 
signal (with n > 2) and the dominant SM backgrounds. We have shown the numbers after the basic 
selection criteria (but before applying any additional cuts) as well as after applying the following 
additional cuts: (i) pf jcts > 50 GeV with p£ ading jct > 100 GeV, and (ii) $ T > 300 GeV. 
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Figure 4: The Jfix distribution for the SISM case with sneutrino LSP (solid lines) and the cMSSM 
case with neutralino LSP (dashed lines) with a similar squark-gluino spectrum. The SM background 
is also shown (shaded region). 



LSP can act as a light DM candidate of mass around 50 GeV while satisfying all the 
existing collider, cosmological as well as low-energy constraints. We also suggest that such 



a scenario can be distinguished from one based on the usual mSUGRA scenario with a 
neutralino LSP, through a study of the same-sign dilepton signals at the LHC, and also 
from the $t spectra in the two cases. We might also be able to put useful bounds on 
the Dirac Yukawa coupling in such scenarios from the invisible decay width of the lightest 
neutral Higgs boson if this gets confirmed with more data at the LHC in near future. 

Acknowledgments 

We thank Florian Staub for useful correspondence. BD would like to thank Haipeng An 
and Rabindra Mohapatra for helpful discussions and collaboration at an earlier stage. BD 
and BM acknowledge the hospitality of the Indian Association for the Cultivation of 
Science (IACS), Kolkata, where a major part of this work was carried out. BD also 
acknowledges the hospitality of the Harish-Chandra Research Institute (HRI), Allahabad, 
during the final stages of this work. This work of BD was supported in part by the 
Lancaster-Manchester-Sheffield Consortium for Fundamental Physics under STFC grant 
ST/J000418/1. BM's work was partially supported by funding available from the 
Department of Atomic Energy, Government of India, for the Regional Centre for 
Accelerator-based Particle Physics, HRI, Allahabad. SM wishes to thank the Department 
of Science and Technology, Government of India, for a Senior Research Fellowship. SR 
acknowledges the hospitality of the Physics Department and the Cluster of Excellence for 
Fundamental Physics 'Origin and Structure of the Universe' at the Technical University 
in Munich (TUM) where a part of this work was done. 

References 

[1] For a review, see e.g., G. Bertone, D. Hooper and J. Silk, Phys. Rept. 405, 279 (2005) 
[hep-ph/0404175]; G. Bertone, Particle Dark Matter, Cambridge University Press (2010). 

[2] J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012). 

[3] M. Drees and G. Gerbier, in Rcf. @. 

[4] For a review, see e.g., G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 
(1996) [hep-ph/9506380]; K. Griest and M. Kamionkowski, Phys. Rept. 333, 167 (2000). 

[5] T. Hebbeker, Phys. Lett. B470, 259 (1999) [arXiv:hep-ph/9910326]. 

[6] T. Falk, K. A. Olive and M. Srednicki, Phys. Lett. B 339, 248 (1994) [hep-ph/9409270]. 

[7] For the latest exclusion limits, see E. Aprile et al. [XENON100 Collaboration], 
arXiv:1207.5988 [astro-ph.CO]. 

[8] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 

[arXiv:1001.4538 [astro-ph.CO]]; D. Larson et al, Astrophys. J. Suppl. 192, 16 (2011) 
[arXiv:1001.4635 [astro-ph.CO]]. 

[9] J. S. Hagelin, G L. Kane and S. Raby, Nucl. Phys. B241, 638 (1984); L. E. Ibanez, Phys. 
Lett. B137, 160 (1984). 

[10] The LEP Electroweak Working Group, http : //lepewwg . web . cern . ch/LEPEWWG. 

[11] R. Bernabei et al. [DAMA Collaboration], Eur. Phys. J. C56, 333 (2008) [arXiv:0804.2741 
[astro-ph]]; ibid, Eur. Phys. J. C67, 39 (2010) [arXiv:1002.1028 [astro-ph.GA]]. 



- 18 - 



[12] C. E. Aalseth et al. [CoGeNT Collaboration], Phys. Rev. Lett. 106, 131301 (2011) 

[arXiv:1002.4703 [astro-ph.CO]]; ibid, Phys. Rev. Lett. 107, 141301 (2011) [arXiv:1106.0650 
[astro-ph.CO]]. 

[13] G. Angloher et al. [CRESST-II Collaboration], Eur. Phys. J. C72, 1971 (2012) 
[arXiv:1109.0702 [astro-ph.CO]]. 

[14] E. Aprilc et al. [XENON100 Collaboration], Phys. Rev. Lett. 107, 131302 (2011) 

[arXiv:1104.2549 [astro-ph.CO]]; J. Angle et al. [XENON10 Collaboration], Phys. Rev. Lett. 
107, 051301 (2011) [arXiv: 1104.3088 [astro-ph.CO]]. 

[15] Z. Ahmed et al. [CDMS-II Collaboration], Science 327, 1619 (2010) [arXiv:0912.3592 

[astro-ph.CO]]; ibid., Phys. Rev. Lett. 106, 131302 (2011) [arXiv: 101 1.2482 [astro-ph.CO]]; 
ibid., arXiv: 1203. 1309 [astro-ph.CO]. 

[16] Z. Ahmed et al. [CDMS and EDELWEISS Collaborations], Phys. Rev. D 84, 011102 (2011) 
[arXiv: 1 105 . 3377 [astro-ph.CO]] . 

[17] E. Armengaud et al. [EDELWEISS Collaboration], arXiv:1207.1815 [astro-ph.CO]. 

[18] V. N. Lebedenko et al. [ZEPLIN-III Collaboration], Phys. Rev. D 80, 052010 (2009) 

[arXiv:0812.1150 [astro-ph]]; D. Y. Akimov et al. [ZEPLIN-III Collaboration], Phys. Lett. B 
709, 14 (2012) [arXiv: 11 10.4769 [astro-ph.CO]]. 

[19] M. Felizardo et al. [SIMPLE Collaboration], Phys. Rev. Lett. 108, 201302 (2012) 
[arXiv:1106.3014 [astro-ph.CO]]. 

[20] P. J. Fox, J. Liu and N. Werner, Phys. Rev. D 83, 103514 (2011) [arXiv:1011.1915 [hep-ph]]; 
J. Kopp, T. Schwetz and J. Zupan, JCAP 1203, 001 (2012) [arXiv:1110.2721 [hep-ph]]; 
C. Kelso, D. Hooper and M. R. Buckley, Phys. Rev. D 85, 043515 (2012) [arXiv:1110.5338 
[astro-ph.CO]]; M. T. Frandsen, F. Kahlhoefer, C. McCabe, S. Sarkar and 
K. Schmidt-Hoberg, JCAP 1201, 024 (2012) [arXiv:1111.0292 [hep-ph]]; P. Gondolo and 
G. B. Gelmini, arXiv:1202.6359 [hep-ph]; R. Foot, arXiv:1203.2387 [hep-ph]. 

[21] J.-F. Grivaz, in K. Nakamura et al. (Particle Data Group), Journal of Physics G37, 075021 
(2010), p. 1312. 

[22] O. Buchmueller and P. de Jong, in Ref. @. 

[23] For global fits of the cMSSM, see e.g., O. Buchmueller et al, arXiv:1112.3564 [hep-ph]; 
C. Strege, G. Bertone, D. G. Cerdeno, M. Fornasa, R. R. de Austri and R. Trotta, JCAP 
1203, 030 (2012) [arXiv:1112.4192 [hep-ph]]; L. Roszkowski, E. M. Sessolo and Y. -L. S. Tsai, 
arXiv:1202.1503 [hep-ph]; J. Ellis and K. A. Olive, Eur. Phys. J. C 72, 2005 (2012) 
[arXiv:1202.3262 [hep-ph]]; A. Fowlie et al., arXiv:1206.0264 [hep-ph]; S. Akula, P. Nath and 
G. Peim, arXiv:1207.1839 [hep-ph]; C. Beskidt, W. de Boer, D.I. Kazakov, F. Ratnikov, 
arXiv:1207.3185 [hep-ph]. 

[24] D. Hooper, T. Plehn, Phys. Lett. B562, 18 (2003) [hep-ph/0212226]. 

[25] H. K. Dreiner, S. Heinemeyer, O. Kittel, U. Langenfeld, A. M. Weber, G. Weiglein, Eur. 

Phys. J. C62, 547 (2009) [arXiv:0901.3485 [hep-ph]]; E. Kuflik, A. Pierce, K. M. Zurek, Phys. 
Rev. D81, 111701 (2010) [arXiv: 1003.0682 [hep-ph]]; N. Fornengo, S. Scopel, A. Bottino, 
Phys. Rev. D83, 015001 (2011) [arXiv: 101 1.4743 [hep-ph]]; L. Calibbi, T. Ota, Y. Takanishi, 
JHEP 1107, 013 (2011) [arXiv:1104.1134 [hep-ph]]; ibid., arXiv:1112.0219 [hep-ph]. 

[26] For a review, see e.g., K. Nakamura and S. T. Petcov, in Ref. [§. 



- 19 - 



[27] P. Minkowski, Phys. Lett. B67, 421 (1977); T. Yanagida, in Workshop on Unified Theories, 
KEK Report No. 79-18 (1979), p. 95; S. L. Glashow, in 1979 Cargese Summer Institute on 
Quarks and Leptons, M. Levy et al. (eds.), Plenum Press (1980), p. 687; M. Gcll-Mann, P. 
Ramond and R. Slansky, in Super gravity, D. Freedman et al. (eds.), North Holland (1980), p. 
315; R. N. Mohapatra and G. Scnjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[28] A. Mazumdar, arXiv: 1106.5408 [hep-ph]. 

[29] R. N. Mohapatra, Phys. Rev. Lett. 56, 561 (1986); R. N. Mohapatra and J. W. F. Valle, 
Phys. Rev. D 34, 1642 (1986). 

[30] C. Arina, F. Bazzocchi, N. Fornengo, J. C. Romao, J. W. F. Valle, Phys. Rev. Lett. 101, 
161802 (2008) [arXiv:0806.3225 [hep-ph]]. 

[31] Z. Kang, J. Li, T. Li, T. Liu, J. Yang, arXiv:1102.5644 [hep-ph]. 

[32] S. Khalil, H. Okada, T. Toma, JHEP 1107, 026 (2011) [arXiv:1102.4249 [hep-ph]]. 

[33] S. Khalil, Phys. Rev. D 82, 077702 (2010) [arXiv:1004.0013 [hep-ph]]. 

[34] F.-X. Josse-Michaux, E. Molinaro, Phys. Rev. D 84, 125021 (2011) [arXiv:1108.0482 
[hep-ph]]. 

[35] L. Basso, O. Fischer, J. J. van dcr Bij, arXiv:1207.3250 [hep-ph]. 

[36] H. An, P. S. Bhupal Dev, Y. Cai and R. N. Mohapatra, Phys. Rev. Lett. 108, 081806 (2012) 
[arXiv:1110.1366 [hep-ph]]. 

[37] P. S. Bhupal Dev, R. N. Mohapatra, Phys. Rev. D 81, 013001 (2010) [arXiv:0910.3924 
[hep-ph]]; Phys. Rev. D. 82, 035014 (2010) [arXiv:1003.6102 [hep-ph]]. 

[38] The CMS Collaboration, CMS-PAS-HIG- 12-020 (2012). 

[39] The ATLAS Collaboration, ATLAS-CONF-2012-093 (2012). 

[40] H. Bacr, V. Barger and A. Mustafayev, Phys. Rev. D 85, 075010 (2012) [arXiv:1112.3017 
[hep-ph]]; ibid., JHEP 1205, 091 (2012) [arXiv:1202.4038 [hep-ph]]; J. L. Feng, 
K. T. Matchev and D. Sanford, Phys. Rev. D 85, 075007 (2012) [arXiv:1112.3021 [hep-ph]]; 
A. Arbey, M. Battaglia, A. Djouadi, F. Mahmoudi and J. Quevillon, Phys. Lett. B 708, 162 
(2012) [arXiv:1112.3028 [hep-ph]]; S. Akula, B. Altunkaynak, D. Feldman, P. Nath and 
G. Peim, Phys. Rev. D 85, 075001 (2012) [arXiv:1112.3645 [hep-ph]]; L. Aparicio, 
D. G. Cerdcno and L. E. Ibanez, JHEP 1204, 126 (2012) [arXiv:1202.0822 [hep-ph]]; 
M. Kadastik, K. Kannikc, A. Racioppi and M. Raidal, JHEP 1205, 061 (2012) 
[arXiv:1112.3647 [hep-ph]]. J. Cao, Z. Heng, D. Li and J. M. Yang, Phys. Lett. B 710, 665 
(2012) [arXiv:1112.4391 [hep-ph]]; N. Karagiannakis, G. Lazarides and C. Pallis, 
arXiv:1201.2111 [hep-ph]; J. Cao, Z. Heng, J. M. Yang and J. Zhu, arXiv:1207.3698 [hep-ph]. 

[41] M. Hirsch, T. Kernreiter, J. C. Romao and A. Villanova del Moral, JHEP 1001, 103 (2010) 
[arXiv:0910.2435 [hep-ph]]. 

[42] E. Ma, Phys. Rev. D 80, 013013 (2009) [arXiv:0904.4450 [hep-ph]]; F. Bazzocchi, 
D. G. Cerdeno, C. Munoz and J. W. F. Valle, Phys. Rev. D 81, 051701 (2010) 
[arXiv:0907.1262 [hep-ph]]. 

[43] S. C. Park, K. Wang and T. T. Yanagida, Phys. Lett. B 685, 309 (2010) [arXiv:0909.2937 
[hep-ph]]; C. S. Fong, R. N. Mohapatra, I. Sung, Phys. Lett. B704, 171-178 (2011) 
[arXiv: 1107.4086 [hep-ph]]. 



- 20 - 



[44] S. Bethke, Eur. Phys. J. C 64, 689 (2009) [arXiv:0908.1135 [hep-ph]]. 

[45] The Tevatron Electroweak Working Group and CDF and DO Collaborations, arXiv: 1107.5255 
[hep-ex] . 

[46] A. Geringer-Sameth and S. M. Koushiappas, Phys. Rev. Lett. 107, 241303 (2011) 
[arXiv:1108.2914 [astro-ph.CO]]. 

[47] M. Ackcrmann et al. [Fcrmi-LAT Collaboration], Phys. Rev. Lett. 107, 241302 (2011) 
[arXiv:1108.3546 [astro-ph.HE]]. 

[48] G. Venanzoni, J. Phys. Conf. Scr. 349, 012008 (2012) [arXiv:1203.1501 [hep-ex]]. 

[49] T. Aoyama, M. Hayakawa, T. Kinoshita and M. Nio, arXiv:1205.5368 [hep-ph]. 

[50] J. P. Lees et al. [BaBar Collaboration], arXiv:1207.2690 [hep-ex]. 

[51] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 108, 231801 (2012) [arXiv:1203.4493 
[hep-ex]]. 

[52] A. Abada, C. Biggio, F. Bonnet, M. B. Gavela and T. Hambyc, JHEP 0712, 061 (2007) 
[arXiv:0707.4058 [hep-ph]]; S. Antusch, J. P. Baumann and E. Fernandez-Martinez, Nucl. 
Phys. B 810, 369 (2009) [arXiv:0807.1003 [hep-ph]]. 

[53] N. Sato et al. (Kamiokande Collaboration), Phys. Rev. D 44, 2220 (1991). 

[54] https : //twiki . cern . ch/twiki/bin/view/CMSPublic/PhysicsResultsSUS 

[55] https : //twiki . cern. ch/twiki/bin/view/AtlasPublic/SupersymmetryPublicResults 

[56] S. Chatrchyan et al. [CMS Collaboration], arXiv:1207.1898 [hep-ex]. 

[57] U. Chattopadhyay and P. Nath, Phys. Rev. Lett. 86, 5854 (2001) [hep-ph/0102157]; ibid., 
Phys. Rev. D 66, 093001 (2002) [hep-ph/0208012]; W. dc Boer, M. Hubcr, C. Sander and 
D. I. Kazakov, Phys. Lett. B 515, 283 (2001). 

[58] See e.g., S. Martin, In "Perspectives on supersymmetry II", G. L. Kane (ed.), p. 1-153 
[hep-ph/9709356]. 

[59] D. V. Forero, M. Tortola and J. W. F. Valle, arXiv:1205.4018 [hep-ph]. 

[60] F. Staub, arXiv:0806.0538 [hep-ph]; Comput. Phys. Commun. 181, 1077 (2010) 

[arXiv:0909.2863 [hep-ph]]; Comput. Phys. Commun. 182, 808 (2011) [arXiv:1002.0840 
[hep-ph]]. 

[61] W. Porod, Comput. Phys. Commun. 153, 275 (2003) [hep-ph/0301101]; W. Porod and 
F. Staub, arXiv: 1104. 1573 [hep-ph]. 

[62] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 176, 367 
(2007) [hep-ph/0607059]; Comput. Phys. Commun. 180, 747 (2009) [arXiv:0803.2360 
[hep-ph]]. 

[63] L. J. Hall, D. Pinner and J. T. Ruderman, JHEP 1204, 131 (2012) [arXiv:1112.2703 
[hep-ph]]; H. Baer, V. Bargcr, P. Huang, A. Mustafaycv and X. Tata, arXiv:1207.3343 
[hep-ph] . 

[64] S. Heinemeyer, D. Stockinger and G. Weiglein, Nucl. Phys. B 690, 62 (2004) 
[hep-ph/0312264]. 

[65] M. Misiak et al, Phys. Rev. Lett. 98, 022002 (2007) [hep-ph/0609232]. 



- 21 - 



[66] S. Bertolini, F. Borzumati, A. Masiero and G. Ridolfi, Nucl. Phys. B 353, 591 (1991); P. 

Nath and R. L. Arnowitt, Phys. Lett. B 336, 395 (1994); F. Borzumati, M. Drees and M. M. 
Nojiri, Phys. Rev. D 51, 341 (1995); H. Bacr, M. Brhlik, D. Castano and X. Tata, Phys. 
Rev. D58, 015007 (1998) [hep-ph/9712305]. 

[67] A. J. Buras, Acta Phys. Polon. B 41, 2487 (2010) [arXiv:1012.1447 [hep-ph]]. 

[68] S. R. Choudhury and N. Gaur, Phys. Lett. B 451, 86 (1999) [hep-ph/9810307]; K. S. Babu 
and C. F. Kolda, Phys. Rev. Lett. 84, 228 (2000) [hep-ph/9909476]; J. K. Mizukoshi, X. Tata 
and Y. Wang, Phys. Rev. D 66, 115003 (2002) [hep-ph/0208078]. 

[69] F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvcstrini, Nucl. Phys. B 477, 321 (1996) 
[hep-ph/9604387]. 

[70] A. Ilakovac and A. Pilaftsis, Nucl. Phys. B 437, 491 (1995) [hep-ph/9403398]]; F. Deppisch 
and J. W. F. Valle, Phys. Rev. D 72, 036001 (2005) [hep-ph/0406040]; ibid., Nucl. Phys. B 
752, 80 (2006) [hep-ph/0512360]; A. Ibarra, E. Molinaro and S. T. Pctcov, Phys. Rev. D 84, 
013005 (2011) [arXiv:1103.6217 [hep-ph]]; W. Abdallah, A. Awad, S. Khalil and H. Okada, 
arXiv:1105.1047 [hep-ph]; A. Abada, D. Das and C. Weiland, JHEP 1203, 100 (2012) 
[arXiv:1111.5836 [hep-ph]]; A. Abada, D. Das, A. Vicente and C. Weiland, arXiv:1206.6497 
[hep-ph] . 

[71] M. Malinsky, T. Ohlsson and H. Zhang, Phys. Rev. D 79, 073009 (2009) [arXiv:0903.1961 
[hep-ph]]; M. Malinsky, T. Ohlsson, Z. -z. Xing and H. Zhang, Phys. Lett. B 679, 242 (2009) 
[arXiv:0905.2889 [hep-ph]]; R. Lai Awasthi and M. K. Parida, arXiv:1112.1826 [hep-ph]. 

[72] K. Ghosh, B. Mukhopadhyaya and U. Sarkar, Phys. Rev. D 84, 015017 (2011) 

[arXiv:1105.5837 [hep-ph]]; Y. Bai, P. Draper and J. Shclton, arXiv:1112.4496 [hep-ph]. 

[73] P. P. Giardino, K. Kannikc, M. Raidal and A. Strumia, JHEP 1206, 117 (2012) 

[arXiv:1203.4254 [hep-ph]]; J. R. Espinosa, C. Grojcan, M. Muhlleitner and M. Trott, 
arXiv:1205.6790 [hep-ph]; ibid., arXiv:1207.1717 [hep-ph]; S. Banerjee, S. Mukhopadhyay and 
B. Mukhopadhyaya, arXiv:1207.3588 [hep-ph]. 

[74] Y. Mambrini, Phys. Rev. D 84, 115017 (2011) [arXiv:1108.0671 [hep-ph]]; ibid., 
arXiv:1112.0011 [hep-ph]. 

[75] P. S. Bhupal Dev, Roberto Franccschini, R. N. Mohapatra, arXiv:1207.2756 [hep-ph]. 

[76] See e.g., H. Baer and X. Tata, arXiv:0805.1905 [hep-ph]. 

[77] G. Bclanger, S. Kraml and A. Lcssa, JHEP 1107, 083 (2011) [arXiv:1105.4878 [hep-ph]]. 

[78] F. del Aguila and J. A. Aguilar-Saavedra, Nucl. Phys. B 813, 22 (2009) [arXiv:0808.2468 
[hep-ph]]. 

[79] F. del Aguila, J. A. Aguilar-Saavedra and J. de Bias, Acta Phys. Polon. B 40, 2901 (2009) 
[arXiv:0910.2720 [hep-ph]]; C. -Y. Chen and P. S. Bhupal Dev, Phys. Rev. D 85, 093018 
(2012) [arXiv:1112.6419 [hep-ph]]; A. Das and N. Okada, arXiv:1207.3734 [hep-ph]. 

[80] S. Mondal, S. Biswas, P. Ghosh and S. Roy, JHEP 1205, 134 (2012) [arXiv:1201.1556 
[hep-ph]]. 

[81] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006) [hep-ph/0603175]. 

[82] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000) [hep-ph/9903282]; 
J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, JHEP 
0207, 012 (2002) [hep-ph/0201195]. 



-22- 



[83] A. J. Barr and C. Gwenlan, Phys. Rev. D 80, 074007 (2009) [arXiv:0907.2713 [hep-ph]]. 

[84] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, JHEP 0307, 001 
(2003) [hep-ph/0206293]. 



- 23 - 



